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Abstract 



We study inhomogeneous solutions of a 3+1- dimensional Einstein-Maxwell-scalar theory. 
Our results provide a holographic model of superconductivity in the presence of a charge 
density wave sourced by a modulated chemical potential. We find that below a critical 
temperature T c superconducting stripes develop. We show that they are thermodynamically 
favored over the normal state by computing the grand canonical potential. We investigate the 
dependence of T c on the modulation's wave vector, which characterizes the inhomogeneity. 
We find that it is qualitatively similar to that expected for a weakly coupled BCS theory, 
but we point out a quantitative difference. Finally, we use our solutions to compute the 
conductivity along the direction of the stripes. 



1. Introduction 



Strongly coupled systems are often found at the heart of spectacular phenomena in con- 
densed matter. Most prominently, they are expected to play a key role in the understanding 
of one of the most exciting scientific discoveries of the last thirty years, high temperature 
superconductivity. While some techniques have been devised over the years to model such 
systems, there is a scarcity of simple, tractable models on which to build our intuition. Re- 
cently, there has been a flurry of activity trying to fill in this gap using various developments 
of the gauge/gravity correspondence. 

The main complexity that arises in strongly correlated electron systems, such as the high 
T c cuprates and the most recent iron pnictides, is the presence of "competing orders" PQ: 
while the system may be a superconductor, it is not a homogeneous one and orders appear 
that are related to the breaking of the lattice symmetries, but are seemingly unrelated to 
superconductivity. Does the ubiquitous presence of such orders in strongly correlated super- 
conductors have a deep connection to the very emergence of superconductivity and/or the 
magnitude of its critical temperature T c ? The answer to this question seems to be crucial for 
identifying the mechanism of superconductivity, at least in the cuprates. It is an important 
empirical observation, derived mainly from numerical calculations of Hubbard models [2], 
that inhomogeneity plays a much bigger role in these systems than in standard, weakly cou- 
pled superconductors. In fact, models have been proposed based on the coexistence of homo- 
geneous superconductivity with charge density (CDW) or spin density waves (SDW), as well 
as models where the superconducting order parameter itself is modulated (pair density waves 
or PDW). Signatures of CDW have been reported in a variety of strongly correlated supercon- 
ductors, most notably the hole-doped cuprates Lai. 6 _ ;r Ndo.4Sr x Cu04, Lai.8-xEuo.2Sr x Cu04, 
and La 2 _ x Ba 3; Cu04 pp. For these materials, the order can be consistently interpreted in 
terms of uni-directional SDW and CDW over a wide range of doping. Altogether, there is 
ubiquitous experimental evidence of ordered inhomogeneous structures in various places of 
the doping-temperature phase diagram of the cuprates, especially in the proximity of the 
superconducting phase. 

Regardless of whether strong coupling and inhomogeneity show us the path to the mystery 
of high T c superconductors or other interesting related systems, it is of great value to have 
simple computable models that incorporate these properties. In this paper we provide one 
model of this kind, in the hope that it might help improve our intuition about these key 
ingredients. 

The AdS/CFT correspondence suggests the equivalence of certain gauge theories at con- 
formal fixed points (CFT) with certain string theories on Anti-de-Sitter (AdS) backgrounds. 
There is a dictionary relating observables in the gauge theory to observables in the string 
theory. One of the key relations is between the curvature of the AdS background and the 
inverse of the gauge, (or rather 't Hooft) coupling of the gauge theory. At large 't Hooft 
coupling and large rank of the gauge group, the D-dimensional CFT becomes dual to gravity 
on D+l-dimensional AdS space, and the correspondence provides a tool to compute ob- 



1 



servables in a strongly coupled theory by solving classical equations of motion on a weakly 
curved background. Further developments of this idea allow one to include perturbations 
away from conformality and develop a richer structure on both sides of the duality. The use of 
this correspondence has allowed, in an unprecedented manner, the study and understanding 
of conformally invariant quantum systems in D dimensions. 

Over the past few years, a growing amount of effort has been devoted to applications 
of the AdS/CFT correspondence to strongly correlated condensed matter systems close to 
quantum critical points. A popular model that includes some important ingredients of a 
realistic system is the holographic superconductor [3E] . It consists of a 3+1-dimensional 
Einstein-Maxwell-scalar theory in an AdS black hole background. At low temperatures, a 
phase transition to a superffuid state takes place [5], in which the scalar field develops a 
non-vanishing expectation value, spontaneously breaking the gauge U(l) symmetry. Physi- 
cal properties such as transport coefficients can be studied and contrasted with other known 
systems like Bardeen-Cooper-Schrieffer (BCS) superconductors. The main conceptual ingre- 
dients of the holographic superconductor are strong coupling and proximity to an underlying 
conformal symmetry, which are believed to be crucial features of cuprates as well. On the 
other hand, inhomogeneity is not accounted for, although it certainly plays a role in the real 
systems. 

In this paper, we would like to raise and attempt to answer a few questions: can one 
take holographic superconductors one step closer to real systems by considering inhomoge- 
neous configurations? Can one account for charge or spin density wave (CDW or SDW) 
backgrounds? Are there any differences in the role of inhomogeneity played in the strongly 
coupled regime as opposed to the weakly coupled BCS case? How do physical properties 
such as T c and anisotropic conductivities compare with experimental data? 

We address these questions by studying inhomogeneous solutions of the holographic su- 
perconductor with a modulated chemical potential of wave vector Q. Throughout this paper, 
we work in the probe approximation, neglecting the backreaction of the Maxwell and scalar 
field on the gravity background. The normal phase then consists of a CDW with the same 
modulation as the chemical potential and a vanishing order parameter. As the temperature 
is decreased the system undergoes a phase transition to a superconducting state with a non- 
vanishing modulated order parameter. We find that various classes of solutions are possible. 
For the two simplest ones, the order parameter is real and the 3-dimensional gauge potential 
vanishes. The periodicity of the order parameter can then be taken to be the same as the one 
of the charge density, or it can be half of it. Here we concentrate on the former possibility, 
but we notice that the latter corresponds to a pair density wave (PDW) [1] . We find that the 
CDW solutions are thermodynamically favored over the PDW solutions, but suspect that 
it may be possible to stabilize the PDW solutions in the presence of a magnetic field. We 
leave a careful study of these solutions for the future. Other solutions that correspond to 
less commonly studied configurations such as spiral pair density waves can also be found. 
The dependence of T c on the wave vector of the modulation (figure [4]) is qualitatively similar 
to the one expected for a weakly coupled BCS system [6]. The two analyses have the same 
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large and small Q asymptotics for T c and both show the presence of a inflection point. The 
main difference is quantitative and concerns the functional dependence of T c on Q in the 
Q — > oo limit. 

We numerically investigate the superconducting state below T c . We find that supercon- 
ducting stripes develop and higher harmonics of the CDW are excited. We verify that this 
state has a lower grand canonical potential (the equivalent of the free energy, but in the 
grand canonical ensemble) than the normal CDW, as expected. Finally, we show results for 
the conductivity along the stripes, verifying the superconducting nature of the inhomoge- 
neous phase. The conductivity in the orthogonal direction is substantially more complicated 
to compute and will be discussed in a future publication. 

The outline of this paper is as follows: in the next section, we introduce the gravitational 
system that we investigate, briefly describe its strongly coupled dual theory, and review 
the construction of a homogeneous holographic superconductor. In section [3j we construct 
an inhomogeneous generalization corresponding to a CDW and show that it undergoes a 
phase transition at low temperatures. The resulting inhomogeneous superconducting phase 
is discussed in section |4| We present results for the grand canonical potential and the 
conductivity in sections [5] and [6j respectively. A discussion of our results is given in section [7j 
Finally, in Appendix [Aj we collect various useful forms of the equations of motion. 



2. The System: Einstein-Maxwell-Scalar Theory in 4D 



The gravitational theory we consider is Einstein gravity in four dimensions with a negative 
cosmological constant A = — 3/L 2 , coupled to a U(l) gauge field and a complex scalar field 
* which is charged under this U(l) symmetry. Here L is a length scale representing the AdS 
curvature radius. The action for this system is 



S 



d 4 Xy/—g 



16nG 



N 



(2.1) 



where D a = d a — iqA a and a, b 6 t,x,y,z. The resulting equations of motion are given 
in (A.l), (A. 2) and (A. 3) of appendix |A} By the gauge/gravity duality [7j, this system on 
an asymptotically anti-de Sitter spacetime is believed to be dual to some large N, strongly 
coupled, conformal gauge theory with a global U(l) symmetry in 2+1 dimensions (for reviews 
see 0)Q 

For concreteness, we will work with a potential 



7(|*|) 



1*1 



(2.2) 



corresponding to a scalar mass m 2 = —2/L 2 . This is the same mass as that of a confor- 
mally coupled scalar in AdS,4. Since this mass is above the Breitenlohner-Friedmann bound 

^ince the £7(1) symmetry on the boundary is global, one is strictly speaking describing a superfluid as 
opposed to a superconductor, but one can imagine weakly gauging the symmetry, making a comparison with 
e.g. BCS theory meaningful. 
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m 2 BF L 2 = —9/4, it does not imply an instability. In the context of the homogeneous holo- 
graphic superconductor, other potentials have been studied (see e.g. [H]), and it has been 
found that the results do not qualitatively depend on the choice of the potential (on the other 
hand see also [ID])- We expect this to be true for the striped holographic superconductors 
as well. 



With this potential, the equations of motion (A.1)-(A.3) possess three linearly indepen- 



dent scaling symmetries under which the various quantities transform as 



X -> X ai X 



(2.3) 



with the charges on given by 





x a 


L 


Q 


A a dx a 




ds 2 


Gn 


Oil 


1 

























1 


-1 


1 





2 





a 3 








-1 


1 


1 





-2 



We can ignore the backreaction of the matter on the geometry if the invariant quantities 
such as Gn | \P | 2 are small. As usual, this can be achieved by sending Gn to zero while keeping 
everything else fixed. Equivalently, using the third scaling symmetry, we can send the gauge 
potential and scalar field to zero while keeping qA and qty fixed, corresponding to a limit 
in which the strength of the gauge interactions is infinitely larger than the gravitational 
interactions. We will exclusively work in this limit in this paper and study the Maxwell and 
scalar fields in a fixed background. The superconductors we would like to describe are at 
finite temperature. So we take the background to be a planar Schwarzschild anti-de Sitter 
black hole which corresponds to studying the dual field theory at finite temperature. We 
choose coordinates so that the line element is 



ds 2 



dz 2 

-h(z)dt 2 + —— + dx 2 + dy 2 
h[z) 



with h(z) 



y 3 ' 



(2.4) 



where z = z Q is the position of the horizon of the black hole. The Hawking temperature of 
this black hole is 

3 



T 



Anzo 



(2.5) 



The equations of motion for the gauge and scalar field (A.l) and (A. 2) with the metric 



corresponding to the line element (2.4) are of course still invariant under our scaling sym- 



metries (provided we assign the same charges to z as to the x a ). We choose to use them to 
set zq — 1, L — 1, and q — 1 to simplify the equations but always plot quantities that are 
invariant under these rescalings. 

For later use, we give the equations of motion in these conventions written in components 
of the gauge field and the real and imaginary part of the scalar field in Appendix A. 
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2.1. Review of the Homogeneous Solutions 



A simple set of solutions to (A.l) and (A. 2) can be obtained jl] by assuming that the fields 
do not depend on x, y and t. In this case it is consistent to set A z = A x = A y = and to 
take ^ = ztjj/\/2 to be real, where the particular definition of ip has been chosen for later 
convenience. With this Ansatz, only two real fields are left: ip and A t , which we denote by 
A for simplicity. 



The equations of motion (A. 4) and (A. 4) on the background (|2.4|) with the assumption 
of homogeneity and stationarity take the forr^ 



hA, 



?P 2 A 



-h 2 ip zz + 3z 2 hip z + [hz - A 2 ) if) 



0. 
0. 



(2.6) 
(2.7) 



A solution is determined once we specify four boundary conditions. For each equation, we 
fix the boundary condition at z — 1, i.e. the horizon of the black hole so that the solution 
is regular. This implies A(l) = and ip z {l) = — ip(l)/3. We choose the third boundary 
condition by requiring that ip has a normalizable profile, which implies ip(0) = 0. The last 
boundary condition determines which ensemble we want to use. To see this, consider a 
solution of the equation of motion ( 2.9[ ) for A and expand it near the AdS boundary as 



A(z) = A (0) + zA w . 



(21 



According to the AdS / CFT dictionary, a gauge symmetry in the bulk theory corresponds to 
a global symmetry in the boundary gauge theory. A^ is mapped into the expectation value 
of the time-component of the current of the global symmetry, while A^ corresponds to a 
source in the boundary Lagrangian. Since A^ can be interpreted as a charge density, its 
canonical conjugate A^ plays the role of a chemical potential. Fixing the charge density or 
the chemical potential corresponds to working in the canonical or grand canonical ensemble, 
respectively. We decide to fix the chemical potential A^ and therefore work in the grand 
canonical ensemble. Notice that in any case, specifying a non-vanishing boundary condition 
for A breaks explicitly the conformal invariance of the boundary theory (which is also broken 
by temperature). 

We first look for solutions in which %p = 0, corresponding to a vanishing order parameter. 
These are dual to the normal, non-superconducting state of the boundary theory. The equa- 
tion of motion for ip is trivially satisfied. We are left with only the simple linear differential 
equation 



A ? 



0. 



(2.9) 



If we choose a chemical potential fi by imposing the boundary conditions ^4(0) = fi we find 
the solution 



A(z) =//(! -z) . 



(2.10) 



2 Once again, we used the symmetries summarized in table [2] to set q = L = Zq = 1. 
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For low chemical potential, or equivalently high temperature, only this normal solution with 
■0 = exists. Beyond a critical value /i c , a new solution appears [oj that has non-vanishing 
condensate, i.e. ip 7^ 0. It has a lower grand canonical potential than the normal solution. 
This solution is typically studied numerically, although some analytical approaches have 
been considered in [TT|12|13IIT^] . The study of the physical properties of the ip 7^ solution 
[3P] reveals that its dual state is some kind of superffuid. The results for the conductivity 
are given in figure [9] and will be thoroughly discussed in section |6j 



3. Inhomogeneous Solutions: the Normal State and T c 



In this section, we begin our study of static inhomogeneous solutions of the system (2.1). 
For previous discussions of inhomogeneous holographic constructions see p2] and p2] , which 
appeared as this work was being finalized. Since we are interested in (static) uni-directional 
charge density waves, we study solutions with non-trivial dependence in the x-direction that 



are independent of y or t. It is then consistent to look for solutions of (A.4)-(A.7) with 



A z = A y = A x = Im\l/ = 0, and the equations of motion with this Ansatz take the fornj^] 

hA zz - ip 2 A + A xx = 0, (3.1) 
-h 2 i) zz + 3z 2 hip z - hip xx + (hz-A 2 )ip = , (3.2) 

where once again 

Similar to what was found in the homogeneous case, we find that for a monochromatic 
chemical potential a phase transition occurs at some critical temperature. In the gauge the- 
ory, we interpret this as a phase transition between a normal and a superfluid/ superconducting 
state in the presence of a charge density wave (CDW). In this section, we study the system 
in the normal phase in which ip = 0, as well as just below T c where the condensate is small. 
This allows us to work at linear order in ip. Since ip enters the equation for A quadratically, 
the equation for the gauge field A in this approximation becomes independent of the scalar, 
simplifying the problem. This allows us to present analytical results for Q = in subsection 



3.2.1 Beyond recovering the homogeneous system, these analytic results provide us with 



some intuition for the effect of inhomogeneities. Since were unable to find analytic results 



for the most interesting range of Q, we provide numerical results for T c in subsection |3. 2. 2 
Working at linear order in the scalar field allows us to find the critical temperature T c at 
which the normal state develops an instability towards the superconducting state tp 7^ 0, but 
is of course insufficient to study the system away from this critical temperature. We return 
to the study of the superconducting solution away from T c in the next section. 



3.1. The CDW Background 



A simple solution of (3.1) and (3.2) is found if ip = and corresponds to the normal non- 
superconducting state. Equation (3.2) is trivially satisfied. For (3.1) we have to specify two 



3 We again set q = L = zq = 1 using the symmetries of the system. 
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boundary conditions. As before, one boundary condition is forced upon us, A(l,x) = 0, if 
we require regularity at the horizon of the black hole. The second boundary condition is the 
physically relevant one since it amounts to specifying the chemical potential (or the charge 
density if we decided to work in the canonical ensemble) of the boundary theory. We decide 
to consider a constant plus a sinusoidal (monochromatic) term imposing 

A(z = 0,x) = /i[(l-6) + 5cos(Qx)) , (3.3) 

where Q is the x-frequency of the sinusoidal term and 5 controls the ratio between inhomo- 
geneous and homogeneous amplitude of the pinning field A(0, x). 5 = corresponds to the 
homogeneous holographic superconductor reviewed in section 2.1 Since (3.1) is linear, we 
can superimpose solutions and look for a solution of the form 

A(z, x) = A (z) + A x (z) cos(Qx), . (3.4) 

The homogeneous contribution to A is then simply given by 

A (z) = - 6)(1 - z) , (3.5) 

while the ^-dependence of A\(z) is governed by 

(l-z^A'l = Q 2 A 1 , (3.6) 

and the boundary conditions are 

A 1 (l,x) = and A x (0, x) = /j,5cos(Qx) . (3.7) 



Although we were not able to find an exact solution of (3.6), we notice that the numerical 
solution is very well approximated by the solution of the simpler equation A'[ = Q 2 Ai, i.e. 



AAz) ~ Ax. 



approx 



z) = fi5 [cosh(Qz) — coth(Q) sinh(Q-2)] 



(3. 



The behavior of the exact numerical solution of (3.6) and the approximate solution ( |3.8 ) 
for various values of Q is given in figure [T] The agreement is excellent both for Q 1 and 
Q > 1, and still reasonably good for Q ~ 0(1). To summarize, the inhomogeneous electric 
background in the normal state (ip = 0) is given by 



Aiormal(£,2:) = A (z) + COs(Qx) A x (z) 

= - 5)(1 - z) + cos(Qx)Ai(z) 



(3.9) 



where Ai(z) is the exact solution of (3.6). 
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Figure 1: The plot shows the numerical exact (continuous red line) and analytically approx- 



imate (dashed black line) solutions of (3.6) for various values of Q. 



3.2. The Instability and the Critical Temperature 

As in the homogeneous case [5], the normal state with ip = becomes unstable for large 
values of fi . To find the value of the chemical potential /i c at which the instability develops, 
we look for marginally stable solutions of (3.2) in the background of (3.9). Once again, 



this is consistent because we are working at linear order in the scalar field ip in the present 
section so that the effect of the condensate on the gauge field is negligible. For the scalar field 
equation of motion, we find it convenient to work in Fourier space and look for a solution of 
the form 



ip(x,z) = '^^ipn(z) cos(nQx) 



(3.10) 



n=0 



The partial differential equation for the scalar field then turns into a system of coupled, 
linear, second order ordinary differential equations for the Fourier coefficients ij) n {z). The ifj n 



satisfy the system (A. 10), which can be conveniently written as 



- V" + (zh(z) -A$)V- A?A U * - 2A A l A 01 + h(z)Q 2 QV = 



(3.H) 
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with \& = (if) Q ,i()i,if)2, ■ ■ ■ ), and the matrices An, A i and Q given by 
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(3.12) 



and 
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(3.13) 



We will now solve the system (3.11) analytically in the limit Q — > 0. Then, in subsection 



3.2.2, we will present the general numerical results and discuss the critical temperature T c . 



3.2.1. Analytical Solution for Q = 

In the limit Q — > 0, the profile of the A field becomes 

lim A(z,x) = A (z) + AAz) cos(Qx) = fi(l - z) [(1 - 5) + 5cos(Qx)} 

Q-s-0 



(3.14) 



where we have used the approximate analytical solution (3.9), which becomes exact in the 



Q — > limit, but have kept the cos(Qx) because it couples different Fourier modes of the 



scalar field. Using (3.14), the system of equations for the Fourier coefficients is 



- + [^K^j I* - ^ 2 (! " z ? [(1 " S? 1 + <* 2 An + 2(5 - 5 2 )A 01 ] ^ = . (3.15) 

Notice that in this limit all the Fourier components have the same ^-dependence, so that it 
becomes natural to look for a solution of the form 

V(z) =v</> A (z) (3.16) 

where v is an eigenvector of the matrix (1 — 5) 2 I + 5 2 An + 2(5 — 5 2 )Aqi satisfying 

[((1 - 5) 2 I + 5 2 A n + 2(5 - 5 2 )A 0l ] v = Av . (3.17) 

The function ip\ then satisfies 

- + (^h(z)) ij x - /i 2 (l - z) 2 \ij x = . (3.18) 



We see that the solution with the largest eigenvalue is the first to become tachyonic. To find 
this eigenvalue, notice that the matrix (1 — 5) 2 I + 5 2 Au + 2(5 — <5 2 )A i is column-stochastic, 
i.e. the entries in each of its columns add up to one. This implies that its largest eigenvalue 
is A = 1. The corresponding eigenvector is 



1 
1 
1 
1 

V ; J 

These are the Fourier coefficients of tt5(x), i.e. for any periodic function 

/(*) 



Oo 

2 



b n 



■ a n cos(nx) + b n sin(nx) 

n=l 

f(x) cos(nx)dx , n > , 

I! 

f(x) sm(nx)dx , n > 1 , 

o 



one finds 



2tt 



dxf(x)S(x)=f(0) 



for 



^ oo 



cos nx 



n=l 



The function ip\=i corresponding to (3.19) satisfies 



- #u + 



Mz) - (1 - z) V 



L^o 



0. 



(3.19) 



(3.20) 

(3.21) 

(3.22) 
(3.23) 

(3.24) 
(3.25) 



(3.26) 



This is the same equation one finds in the homogeneous case. So we see that in the limit 
Q — > 0, independent of 5, we find the same critical temperature as in the homogeneous 
system. We will call this temperature T c ^. This property can be understood intuitively as 
follows. For small Q the x derivatives in (3.2) are irrelevant, so the system is effectively 



equivalent to a collection of homogeneous systems each one at a fixed x, with chemical 
potentials ranging from fi(l — 25) for x = tt/Q, up to fi for x = 0. Therefore, for any 5 there 
is an (isolated) homogeneous system with chemical potential /i , sitting at x = 0, which will 
go through a phase transition precisely at {T/qn) c h . In addition to this result for T c , which 
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Figure 2: We show the z profiles of ifj n for the first few n and two different values of Q at the 
phase transition (T/(q/i) = 0.0269 for Q = 4 and T/(q/j) = 0.0575 for Q = 1/10). Dashed 
lines with longer dashes correspond to larger n. It is evident that ip n are suppressed for 
larger n and more so for larger Q. This results are obtained solving the (truncated) linear 



system (3.11), therefore the normalization of ip n is arbitrary. 



was expected, we can learn a few more things from this analytical computation. For non-zero 
Q, the higher Fourier modes are very costly and will be suppressed. One thus sees that the 
effect of the coupling to the modulated gauge potential tends to excite all the Fourier modes 
democratically while the gradients tend to suppress the high-n modes. This suggests that it 
will be a good approximation to set to zero the Fourier modes above some n max which will 
depend on Q and the desired degree of accuracy. 



3.2.2. Numerical Results for T c 



In the following we present the numerical solutions of the system (3.11) and discuss the 
critical temperature. Let us start by describing the numerical method we use. As we 
noticed, the effect of Q in (3.11 ) is to suppress ip n for large n. For a given momentum Q and 
precision e, there is an n max (Q,e) such that all ip n for n > n max (Q,e) contribute to if) less 
than e and can be neglected. This is evident in figure [2] where the first ip n are plotted for 



two different values of Q. This means we are allowed to truncate the hierarchy (3.11), and 
we are left with a finite number of n max (Q, e) coupled linear ordinary differential equations, 
which we solve numerically in Mathematica. In all the numerical results we present we have 
fixed the precision to be 10 -3 . Then, for example taking Q = 1/10, we need to keep as many 
as 14 modes, while for Q = 4 only the first 3 are necessary to achieve the required precision. 



Once we have solutions of (3.11 ), the most interesting physical property we can compute 



is the critical temperature T c at which the normal state becomes unstable. Let us pause to 



explain how we extract this quantity from the solutions of (3.11). Let us choose some fixed 
< 5 < 1 and Q > 0. Then our theory has four parameters (L,q, zq, fi). Using the three 
symmetries in table |2j we can set three of them to an arbitrary value and we once again set 
L = zq = q = 1. Then we numerically solve (3.11) by integrating from the horizon at z = 1 
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Figure 3: The plots show the z profiles of ip n and A n obtained from the (truncated) full 
system (A. 14), for the first few n and 5 = 3/5, Q — 1/4 and T/(g/i) = 0.0556 which is 
very close to the critical temperature. Again, dashed lines with longer dashes correspond to 
larger n. The A n profiles show that close to T c only the zeroth and first harmonics of A are 
excited, which confirms that the backreaction of ip is negligible. 
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Figure 4: The plot shows the critical temperature as function of Q for 5 = 0.2, 0.4, 0.6, 0. 
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out to the boundary at z = 0. To set the initial conditions for this integration, we solve 
the equations as a power series near the horizon. If we impose regularity at the horizon, 
we will have one free coefficient for Aq, one for A\, and one undetermined coefficient for 
each of the n max + 1 Fourier modes of ip giving a total of n max + 3 coefficients. Since the 
equation for if) is linear, the overall normalization is unphysical and we can set the coefficient 
for, say, i[)q to unity. The remaining n max + 2 coefficients are chosen so that ^„(0) = for 
all rimax + 1 modes and SA (0) = (1 — S)Ai(0). We then have n max + 2 equations for the 
n max + 2 coefficients, which fixes all of them. In practice, we do not integrate to z = but 



to z — 10 and match the solution there to a power series of the solution of our (3.11 ) near 
the boundary to extract A o (0), A\(0) and V'n(O). The resulting \i = y4 (0)/(l — 5) is then /i c . 
The solution is unique if we require our Fourier coefficients ip n (z) to have no zeros between 
the horizon and the boundary. The solutions with one or more zeros have higher [i and thus 
lower T and are not relevant for us. This fi is of course not invariant under the symmetries 
in table |5J but the combination T/(g/i) is, so we can compute it for an arbitrary choice of 
those rescalings. This quantity is useful since variations of it can be interpreted, upon an 
appropriate rescaling, as variations of T at constant /i and q. This is the quantity we will 
use in all of our plots involving temperature. 

In figure [4] we have plotted the critical temperature, obtained as explained above, as a 
function of Q for a few different values of 5. We notice the following features: 

• The limit of (T/g/i) c for Q — > is independent of S and it corresponds to the critical 
value of a homogeneous superconductor with chemical potential fi, which we denote 
by {T/qn) c .h- This nicely agrees with our computation and discussion in the last 
subsection. 



The limit of (T/g/i) c for Q — )• oo depends on 6 simply as 



hm - =(!-<$) - „ =(!-*) _ . (3.27) 



T 



c.h 



This can be explained as follows. For very large Q, Ai(z) goes to zero for any z ^ 0, 
as can be seen in figure [Tj Then the background experienced by ip is that of the x- 
homogeneous field A = /i(l — S)(l — z) almost everywhere. Furthermore, as explained 
before, for Q — > oo only ipo is non-zero. Hence the system is equivalent to a homoge- 
neous one with chemical potential /j,(l — S) and therefore experiences a phase transition 
at (1 - S)(T/qn) Cjh . 

One might ask how (T/qfi) c (Q,5) reaches its asymptotic value for Q — > oo. Since Ai 
vanishes exponentially in Qzq in the Q — > oo limit for any z ^ 0, we expect it to behave 
as 

lim (—)~a-S)(—) + e~ Qzo S ( — ] , (3.28) 



13 



plus terms of higher order in e ® z ° which become negligible for large Q. We confirm 
that this behavior is consistent with that of our numerical results by performing various 
fits. 

It is interesting to compare the functional dependence of T c on Q for our striped holo- 
graphic superconductor with that of a weakly coupled BCS superconductor studied in [6] 
using a mean-field approximation. Since we choose a simple consistent ansatz in which the 
order parameter is real, our analysis does not account for phase fluctuations and phase co- 
herence, so we should compare with the continuous line in figure 1 of [6]. The two results 
look extremely similar. Both the Q = and Q = oo asymptotics are the same and in 
both cases there is an inflection point at finite Q. The main difference consists in how the 
Q = oo limit is reached. Expanding (4) of [6] for large Q and translating into our notation 
for convenience, we find 

BCS : lim T c ~ (1 - 5)T c , h + T-^-—T c h + O ( 1 ) , (3.29) 

logQ \(\ogQ) 2 J 



to be contrasted with the exponential behavior of (3.28 ). It would be interesting to generalize 



our analysis to include phase fluctuations and compare the results with the expectations of 



4. Inhomogeneous Solutions: the Superconducting State 



In this section we study the full system of coupled non-linear partial differential equations 



in (3.2) and (3.1). We start by expanding both A and ip m Fourier modes as 



i>(x,z) = y^^ n (z) cos(nQx) , 

n=0 
oo 

A(x,z) = A n {z) cos(nQx) , 



(4.1) 
(4.2) 



n=0 



and plugging this back into (3.2 ) and (3.1 ). Then, by projecting onto each Fourier mode, we 
obtain a hierarchy of coupled non-linear ordinary differential equations for A n (z) and ip n {z). 
As before, we can truncate the expansion at some n max (Q, e), which we again choose so that 
the modes ipnm ax and A nmax are a factor of 10~ 3 smaller than the ip n or A n with the largest 
amplitude. The hierarchy is too long to be written out here, but for the convenience of the 



reader we include a few lines of Mathematica code (A. 14) in the appendix that can compute 



it. To solve the hierarchy, we again integrate it from the horizon out to the boundary. The 
boundary conditions at the horizon are again obtained from a power series of the solution near 
the horizon. After setting A n (l) = to ensure regularity, one has 2n max + 2 undetermined 
coefficients. We require ip n (0) = for the (n max + 1) Fourier modes of the scalar field in order 
to have a normalizable profile. We also impose A (0) = — S), Ai(0) = fiS and A n (0) = 
for 1 < n < n max . This provides 2n max + 2 equations for the 2n max + 2 Fourier coefficients. 
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Qz =0.25 £=0.4 T/q//=0.03704 Qz =0.25 £=0.4 T/q^=0 .03704 




Qz =2. £=0.4 T/qyu=0.02919 Qz =2. £=0.4 T/q^=0.02919 




Figure 5: The plots show the z, x profile of A and ip for Q = .25 and Q = 2. The striped pat- 
tern of condensation of the bulk scalar field ip is evident in both cases, but more pronounced 
for lower Q's. 
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Figure 6: The left and right plots show the first few ip n and A n (z), respectively, obtained 
from the (truncated) full system (A. 14) for Q = 1/4, T/(q/j) = and 5 = 3/5. Longer dashes 
correspond to larger n. The Fourier modes A n for n > 1 are excited by the backreaction of 
the scalar condensate. 




Figure 7: The dependence of the superconducting order parameter 2 = d z ip(0,x) on tem- 
perature and x is plotted for Q = 1/4 and Q = 2 with 5 = 2/5. We plot two periods in x to 
better visualize the superconducting stripes, which have the same periodicity as the CDW. 
Notice that for small Q the system is superconducting only in the stripes since the order 
parameter is very small in between. For larger values of Q, the system is stiffer in the x 
direction and therefore the stripes are less pronounced. 



16 



Again, we only integrate to z — 10~ 8 and match to a power series expansion of the solution 
near the boundary to extract these quantities. As before, we choose the solution that has 
no zeros in the z-profile, which gives rise to a unique solution. The fact that the truncation 
is a good approximation can be seen in figure [6] where we show a solution of the full system 
(A. 14) for A n (z) and ip n for the first few n. It is clear that both A n and if) n become smaller 
and smaller as n increases. 

In figure [5] we show if)(z, x) and A(z, x) as obtained from our numerical computation with 
the specified values of T, 5 and Q. At some fixed subcritical T < T c , as Q is gradually in- 
creased (keeping everything else fixed), there is a smooth crossover between two qualitatively 
distinct types of profiles of ip. For Q < 1, ijj has a sharp striped condensation pattern in the 
bulk: the scalar field has a large expectation value in correspondence of the maxima of the 
CDW, while it is vanishingly small in between as can be seen in the top half of figure [5| This 
behavior carries over to the order parameter of the boundary theory, i.e. (^(x) = d z ip(0, x), 
as we show on the left of figure [7} The boundary system has sharp superconducting stripes, 
well separated by stripes of normal phase. The situation changes continuously as we increase 
Q. For Q > 1, the system is stiff enough in the x direction that bulk condensation takes 
place everywhere at approximately the same temperature, although modulations are still 
clearly visible (bottom of figure [5]). This translates into the fact that the boundary system 
is well inside the superconducting phase everywhere, but the order parameter is modulated 
as we show on the right of figure [7j 



5. The Grand Canonical Potential 

In this section we study the grand canonical potential Q for the system. This is the Legendre 
transform of the free energy F, which is the thermodynamic potential in the canonical 
ensemble. Before proceeding, let us be pedantic and spell out this difference in detail. The 
natural variables for F are (T, p), while those for Q are (T, p), where p is the chemical 
potential for p. Using the AdS/CFT dictionary, in our computation plays the role of 
the chemical potential p, since it perturbs the Lagrangian of the boundary theory, while 
A^- 1 ' plays the role of the (charge) density p, since it is the expectation value of the time 
component of a current, i.e. a density. We have decided to give our results in terms of T / (gp) 
which is invariant under all the symmetries in table [2} All our plots are directly comparable 
with any other convention for fixing three out of the four parameters {q, Zq, p, L} (in the 
explicit computations we chose q — L — Zq — 1). Variations of T/(gp) can be naturally 
interpreted as variations of T at constant p and q. This means that we are specifying T 
and p, while the entropy and p are dependent variables. In particular notice that if we fix 
p and vary T, in general p will vary as can be seen by any explicit solution of the equation 
of motion for A. To summarize, we are working in the grand canonical ensemble^] and we 

4 Of course one could insist on using T/(q/j,) and to work with the canonical ensemble, i.e. to keep p fixed. 
In that case variation of T/(qp) would be interpreted as some simultaneous variation of T and fi, which has 
not a very transparent physical interpretation. In the canonical ensemble the results should be naturally 
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should use the grand canonical potential Q. 

There are two competing states, the normal state if) = and the superconducting state 
if) ^ 0. We show that below T c the grand canonical potential of the latter is lower, confirming 
that a phase transition takes place. We find that the order parameter varies continuously as 
the temperature crosses T c , while its first derivative does not. Hence the phase transition to 
the superconducting state is of the second order. 

By the AdS/CFT dictionary, Q/kT is the negative of the on-shell action of the Wick- 
rotated solution with Euclidean time compactified on a circle of circumference 1/kT. Very 
often this quantity is infinite and one needs to renormalize it. We present this derivation 
for our inhomogeneous holographic superconductor following [17] . but see also [T8fl9|l20] for 
related computations. We are interested in the action for A ancQ \1/ assuming a static and 
^/-independent solutions. We will assume that \P is real as it is the case for all our solutions. 
This can be straightforwardly derived from (2.1) 

(d x A) 21 



S 



d x- 



1 



(d z A) + 



h 



h{d z ^y + (d x ^y - ^ 



2^2 



A 2 z 
~h 



+ 2 



(5.1) 



To evaluate S on shell it is convenient to integrate by parts and use the equations of motion. 
This leads to 



S = J d 4 x [e.o.m.'s] + S, 



on— shell 



s 



on— shell 



z=0 



4**' - -AA' 
2 2 



z 



~ 1 ^ 



(5.2) 
(5.3) 



where we have used the fact that A(l,x) = h(l) = while *&(l,x) is finite. S^n-sheii has 
divergences coming from the boundary term at z — 0. We therefore introduce a cutoff at 
z = e to regularize it. The on-shell action can then be expanded as 



S, 



on— shell 



d 3 x 



[VW(x)]' 



+ finite terms , 



where we have used the fact that close to the boundary \I/ is given by 

(z, x) = ¥ 1) (x)z + ¥ 2 \x)z 2 + ... . 
This divergence can be subtracted by adding the following boundary counter term 



S c .t = 



d 3 x\& 



Hence the renormalized grand canonical potential is 



n 



dxdy 



2 = 



+ 



A 2 

dxdydz- — r^ 2 
hz 2 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



plotted as function of T/q^Jp. 

5 For this derivation formulae are simpler if we use ^ instead of ip. Once we obtain the final result, we 
will convert it to i\) = V2Re\E'/z. 
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Figure 8: The plot shows the grand canonical potential Q (multiplied by the appropriate 
factors to make it invariant under the symmetries in table [2]) as function of T/(g/z). Below 
the critical temperature both the normal and the superconducting states exist. The latter 
always has a lower Q. 



We can simplify this formula when considering only normalizable profiles of i.e. = 0, 
which is the case for all the solutions considered in this paper. The average grand canonical 
potential per unit volume then becomes 



n 
v 



Q_ 

2tt 



2 = 



dx l A A' + Q- 
2 2vr 



A 2 

dxdz- — -ip 2 
hz 2 



(5.8) 



where we have reverted to the notation ip = \/2Ke^ / z. The expression for Vt further sim- 
plifies in the normal state where ip = 0. There we can compute (5.8) analytically using our 
approximate solution for the normal state profile of A 



A(z,x) = A + Ai cos(Qx) , 
A (z) = fi(l-8)(l-z), 
Ai(z) ~ fi5 [cosh(Qz) — coth(Q) smh(Qz)] 



We find 



n 

V 



"2" 



;i-5) 2 + ^ 2 Qcoth(Q) 



(5.9) 
(5.10) 
(5.11) 



(5.12) 



For the superconducting state we have to resort to numerics. It is straightforward, given 



solutions for A and ip, to numerically evaluate (5.8). In figure 8l we show the temperature 
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dependence of Q/V for several values of S for Q = 2 as an example. In all cases there is a 
critical value of T j (qfi) below which two possible solutions exist, the normal one with ip = 
and the superconducting one ip ^ 0. For all values of Q and 5 we have studied, the latter 
has a lower Q and is therefore favored. Notice that the numerical computation of Q for the 



normal state is well approximated by the analytical formula (5.12). 



6. Conductivity 

In this section, we study the electrical conductivity of the boundary theory in the presence 
of a CDW and both normal and superconducting states. Our starting point are the Maxwell 



equations (A.4)-(A.7). We want to add a small homogeneous electric field on top of the 
solutions we have found in sections [3] and [4] and study the linear response of the system. 
Because our solutions are inhomogeneous in the x-direction, an electric field in the x-direction 
sources five other independent perturbations even at linear order. The computation of a x is 
therefore much more complicated than the computations of conductivity currently performed 
in the literature. We hope to present results for it in a future publication. 

The calculation simplifies for a y . In this case, we want to probe the system with a small 
homogeneous electric field in the ^-direction, which we will realize by switching on a time 
dependent A y , and read out the linear response. All the solutions we have presented are 
homogeneous in the y-direction and there are no d y derivatives to contract the perturbation 
A y with. Therefore, at linear order in A y , we only have to solve the Maxwell equation for 
A y (u, z,x) (which is now complex because of the Fourier transform of time into frequency) 



on the backgrounds computed in the previous sections. Using (A. 7) we find 



hd z {hd z A y ) + hd 2 x A y + (u 2 - ip 2 h) A y = 0. (6.1) 

Notice that, since the background is static, the Fourier modes A y (co, x, z) with different 
frequencies u decouple from each other. The optical conductivity in the ^-direction is given 
by 

where we have expanded A y (x, z) around the AdS boundary z = as 

A y (x,z) = Af(x) + zA y l \x) + ... . (6.3) 



Let us discuss the two boundary conditions for (6.1). In the x-direction, we impose peri- 
odic boundary conditions. In the z-direction one boundary condition is simply the overall 
normalization. As long as A y is small enough we can work at linear order in A y . The nor- 
malization is then irrelevant because it cancels in the definition of a y . The second boundary 
condition has to be imposed at the horizon z = z = 1 and is dictated by the requirement 
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Qz =0.5 5=0.4 T/r c =l Qz =0.5 5=0.4 T/T,=0.85 




Figure 9: The panel shows the y-conductivity a y as function of x and of the frequency oj for 
different values of the temperature and Q = .5 and 5 = .4. Since Q is small, the system is 
easily understood in terms of the homogeneous results for the conductivity. A gap opens in 
correspondence of the superconducting stripes x = + 27m, where the order parameter is 
large (as seen e.g. in figure [71), while a y (ui, x) is constant in between the stripes where ip ~ 0. 
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Qz =4 5=0.6 T/7>1 Qz =4 5=0.6 T/7>0.7 




Figure 10: The panel shows the ^/-conductivity o y as function of x and of the frequency uj 
for different values of the temperature and Q = 4 and S = .4. The profile is qualitatively 
very different from the one for Q < 1 in figure |9j A gap opens up everywhere in x, but the 
large u asymptotic value is reached very differently at different points. 
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Figure 11: The plots show a comparison of the conductivity for the homogeneous holographic 
superconductor (dashed line) with an x = slice of the inhomogeneous conductivity a y (u, 0) 
for two different temperatures T/T c = 0.85,0.65. On the left we show Q = 0.5 while on 
the right Q = 4. Once the frequency u is rescaled by the value of the corresponding order 
parameter \fO~^ the Q = 0.5 inhomogeneous conductivity becomes almost indistinguishable 
from the homogenous one. The Q = 4 case shows a qualitatively different behavior. 



of causal behavior, which implies ingoing boundary conditions [21]. In our case this i^] 
A y (u>, z,x) oc (1 — z) +luJ / 3 + . . . near the black hole horizon z = 1. 



We have numerically solved (6.1) and computed a y using (6.2). The results for different 

Let us comment on these results 



values of Q and T / (g/i) are shown in figures M and 10 



starting from the normal CDW background studied in section [3] Since A is a vector boson 
of an Abelian gauge theory, it does not couple to itself and knows about the background only 
through its coupling to the scalar field. In the normal phase in which the scalar field vanishes, 
the conductivity in the presence of a CDW is then the same as the one in a homogeneous state, 
which is <y(u) = 1, in the limit in which the backreaction on the geometry is neglected [3H] 
(for earlier discussion see [22]). 

Things become more interesting below T c where a condensate has developed. In figures 



9] and 10, we show how a y (ui,x) changes as the temperature is lowered for two different 
values of Q. Again, two qualitatively different behavior appear depending on the value of 
Q. For Q < 1, the system effectively behaves as a collection of independent homogeneous 
holographic superconductors, each one at a fixed value of x. This can be formally seen in 
the analysis of subsection |3 . 2 . 1 1 of the Q — > limit. Therefore the conductivity for Q = .5 in 



6 In order to see this, plug the ansatz A y (ui, z, x) = (1 — z) a ^ A y (x) + . . . into ( |6.1[ ) and expand around 
z = 1. At leading order one finds a = ±iw/3. The right sign for ingoing boundary condition, is found 
by Fourier transforming back to the time domain (this defines our conventions of the sign of the Fourier 
transform) 



A y {t, z,x) — 



— t 
2tt 



■[t±log(l- 



'■*>WA y (x) + .. 



(6.4) 



Since log(l — z) decreases going towards the horizon, the ingoing boundary condition is given by the upper 
sign, i.e. A y (uJ, z, x) oc (1 - z) +%UJ / 3 + 
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figure [9] can be qualitatively understood in terms of the conductivity for the homogeneous 
case. A gap forms for low frequencies in correspondence of the superconducting stripes 
x = + 27m, while in between the stripes, where ip ~ 0, the conductivity is the same as 
that of the normal state, i.e. constant. As Q becomes larger than one (keeping all other 
parameters fixed), there is a smooth crossover to a qualitatively different behavior. We plot 
the optical conductivity for Q = 4 in figure 10 The profile is dramatically different from the 
one in figure |9j A gap opens almost simultaneously everywhere in x due to the increased 
stiffness in that direction. Interesting resonant patters are visible as the large u asymptotic 
value is reached. To make the new features of the Q > 1 conductivity even more evident, 



in figure 11 we show slices of a y (u,x) at the center of the stripe (i.e. x — 0) and compare 



them with the homogenous conductivity for two subcritical temperatures. Upon rescaling 
uj by the corresponding value of the order parameter y/O?,, it is clear (left plot of figure 11) 
that the Q = 0.5 conductivity agrees well with the homogeneous one. On the right plot of 



figure 11, where Q = 4, a y (uj,0) presents extra oscillations which are completely absent in 
the homogeneous case. It would be desirable to better understand the physics behind the 
Q > 1 conductivities. 



7. Conclusions 

In this paper we have studied a holographic model of superconductivity in the presence of a 
CDW. Starting from a monochromatic CDW in the normal state with wave vector Q, after 
the phase transition we find superconducting stripes coexisting with the original CDW plus 
its higher harmonics. The dependence of the temperature of the phase transition T c on Q 
has some similarities with what is expected from a weakly coupled BCS computation [6].T C 
decreases monotonically as Q increases, and asymptotes some constant value. The main 
difference is a steeper functional dependence for Q — > oo for the holographic model. These 
conclusions are derived in the absence of phase fluctuations of the condensate, which is a 
consistent assumption at the level of our equations. In [6J it was argued that once phase 
fluctuations are included, there is an optimal Q for which T c has a maximum. It would 
be very interesting to verify this expectation by generalizing our computation to include 
these fluctuations. A first step would be a stability analysis of our solutions under small 
fluctuations. 

To characterize the superconducting phase, we have presented results for the conduc- 
tivity CT y in the direction parallel to the superconducting stripes. The conductivity in the 
perpendicular direction, a x is more complicated to compute because, even at the level of 
linear response theory, several modes are excited. On the other hand, this is a very inter- 
esting observable for a couple of reasons. First, it tells us about the degree of anisotropy of 
the system, which can be directly compared with experiment. Second, the knowledge of o x 
carries information about the correlation between stripes and the nature of proximity effects 
in a holographic model. 

Our analysis has been performed in the probe limit in which the backreaction of the 
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Maxwell and scalar field on the gravitational background are neglected. Although this is 
expected to capture the relevant physics of the system, it might be interesting to study what 
happens when gravity is dynamical. One might worry that inhomogeneities might trigger a 
gravitational collapse making the system unstable. We have two comments in this regard. 
First, it is well known that stability conditions in AdS are different from and often less 
intuitive than those in flat space. Second, since q 2 L 2 /Gn is still a free parameter, one can 
always make the time-scale of the instability arbitrarily large implying that the system is 
effectively stable. 

There is significant information we could acquire by including the equations for gravity. 
It would for example be interesting to see whether the DC conductivity perpendicular to the 
stripes in the normal state has a delta function since translational invariance is broken (as 
e.g. in [23]). One could also see how the findings of (21] about the absence of a hard gap 
change due to the CDW. 

Our focus in this work was to study how superconductivity interacts with a CDW and 
what the physical consequences of this coexistence are. We did not try to spontaneously 
generate the inhomogeneity of the CDW solution, but chose to source it with a modulated 
chemical potential. One could extend our analysis to a model in which the CDW solutions 
are generated dynamically, e.g. via a spontaneous symmetry breaking along the lines of [16J. 

We would like to conclude mentioning that we have found several instabilities of the 
normal state of the holographic superconductor in the presence of a CDW beyond the one 
studied in this work. Another simple one is towards a superconducting state in which 
the order parameter is modulated with half the frequency of the CDW. This construction 
has been proposed in the condensed matter literature to account for various properties of 
the cuprates and is known as pair density wave (PDW) [lj. A thorough discussion of the 
holographic realization of PDW states will be presented elsewhere. 
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A. The Equations of Motion in Components 



In this appendix, we collect the equations of motion derived from the action (2.1) 
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The scalar equation is [3] 

1 



1 



=D a {^-gg ah D b m) + ~V\M) = 0. 



the equation for the U(l) gauge field is 

1 



-g 



d a {V=gF ab ) = iqg ba - 



(a.i; 



(A.2) 



and Einstein's equations are 

1 3 
Rab — -^9abR — j2^ ab = S^Gn 
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V 



F ac Fb C — -g a bF cd F cd 



- 9ab D a y(D b yy + [D a v(D b vy + a^b\- g ah v{\^\) 



. (A.3) 



Breaking the complex scalar field up into its real and imaginary part as $ = (0 + ix)/y2 
and using the metric (2.4), the equations of motion for the two real scalar fields <f> and \ 
become 



-g 



,d a (V^A a ) + 2A a d a 
! 

d a (>/=£A°) + 2A a d a 



X 







0. 



Analogously, assuming no dependence on the y-coordinate, as it is the case for all solutions 
considered in this paper, the Maxwell equations become 

" 2 [hd 2 z A t + d 2 x A t }-A t (<p 2 + x 2 ) = 

z 2 

Xd t (j> - (j>d t x + jo (d t d x A x + hd t d z A z ) , 

v2 



Z 



(A.4) 



z 
J? 



d^A z — —d^Az 



A z (<P 2 + x 2 ) 



Xd z (j) - <pd z x + ( d z d x A x - -d t d z A 



(A.5) 



z 

1? 



(hd 2 A x + h'd z A x ) - -d 2 A x 



A X (<P 2 + X 2 ) 



Xd x <p - <f)d x x + -jt 2 ( hd z d x A z + h'd x A z 



z 

L2 



d z (hd z A y ) - -d 2 t A y + d 2 x A y 



A„ 



h 

l2 i ,,2 



d t d x A t J , 

+ x 2 ) = o. 



(A.6) 
(A.7) 



In section |3| we studied the equation of motion of if) (3.2) in the normal state background 
given by 



.4 



normal 



(z, x) = A + cos(Qx)Ai(z) = //(l — 8)(1 — z) + cos(Qx)A\{z) 



(A.8) 
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where A± is the solution of (3.6) with the boundary condition Ai(0) = 5fi and A\(l) = 0. 
Upon expanding ip in Fourier modes as 



ip(x,z) = j> n ( z ) cos(nQx) 



n=0 



we find an infinite system of linear coupled ordinary differential equations 



+ 



://(:)- [AI + -AI 
zh(z) -[Al + \A\ 



4 



zh{z) — 



K + \A\ 



i[> - AqA^ - -A\^ 2 = 



^ 1 -A A x {2% + f 2 ) 
A?(V>i + V 3 )+/*QVi = 0, 



--^(2^ + ^4) + /i(2Q) 2 ^2 = 



-< + 



0/1(0) — 



^ + ^? 



- A Ai(V> n -l + Vn+l) 



(A.9) 



--A?(^- 2 + ^+2) + /i(^Q) 2 ^2 = , (A.10) 



where Aq(z) and A\(z) are the only two Fourier modes of A being non-zero. As an aside, 
we notice that in the case A = the odd modes ip2n-i decouple from the even modes ip2n- 
The full hierarchy for both A and ip is obtained by plugging the expansions 



ip{x,z) 
A(x,z) 



y^^njz) cos(nQx) 

n=0 

00 

A n (z) cos(nQx) 



(A.ll) 
(A.12) 



n=0 



into (3.2) and (3.1). The final form is too long to be written down here, but we give below 
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the Mathematica commands that produce any given order, in this case five 

nmax = 5; (A. 13) 

h[z_] = 1 - z 3 ; 

0[z_, x_] = Sum[A[n][z]Cos[nQx], {n, 0, nmax}]; 
V>[z_, x_] = Sum[-^[n] [z]Cos[nQx], {n, 0, nmax}]; 

eqa[Q_] = h[z]D[0[z, x], z, z] + D[</>[z, x], x, x] — ip[z, x] 2 /z 2 0[z, x]/ /TrigReduce; 
eqp[Q_] = -h[z]z 2 D[h[z]/z 2 D[^[z,x],z],z] - 2h[z]/z 2 ^[z, x] 

— h[z]D[0[z, x], x, x] — 0[z, x] 2 ?/>[z, x]/ /TrigReduce; 

Do[{eqpQ[n][Q_] = Collect[lf [n == 0, Q/2/PiIntegrate[eqp[Q], {x, 0, 2Pi/Q}], 

Coef f icient[eqp[Q], Cos[nQx]]], {-(//'[n], i/j'[n], 4>[n}}, Simplify]}, {n, 0, nmax}] 

Do[{eqpQ[n][Q_] = Collect[lf [n == 0, Q/2/PiIntegrate[eqa[Q], {x, 0, 2Pi/Q}], 

Coef f icient[eqa[Q], CosfnQx]]], {A"[n], A'[n], A[n]}, Simplify]}, {n, 0, nmax}] 
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